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METHOD FOR PROVIDING ERROR INFORMATION RELATING TO 
INCONSISTENCIES IN A SYSTEM OF DIFFERENTIAL EQUATIONS 

5 

Cross-Reference to Related Application : 

This application is a continuation of copending International 
Application No. PCT/DE02/01878 , filed May 23, 2002, which 
designated the United States and which was not published in 
10 English. 

Background of the Invention : 
Field of the Invention : 
The invention relates to a method for providing error 
15 information relating to inconsistencies in a system of 
differential equations . 

It is known in the prior art to describe a technical system 
with systems of differential equations. With the assistance of 

2 0 numerical methods using a computer system, or of an analog 
computer, it is possible to simulate the systems by 
calculating solutions to the systems of differential equations 
for specific initial or boundary conditions. Examples for the 
description of technical systems and typical simulation 

25 methods can be taken from standard works of the relevant field 
of specialization. Detailed instructions therefor are 
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specified, for example, in G. Schmidt, "Grundlagen der 
Regelungstechnik" ["Principles of Control Engineering"] , 2nd 
edition, Springer, Berlin, 1987, in Unbehauen, 
"Regelungstechnik I" ["Control Engineering I"], 6th edition, 
5 Friedr. Viehweg & Sohn, Braunschweig/Wiesbaden, 1989, in 

Unbehauen, "Regelungstechnik II" ["Control Engineering II"] , 
5th edition, Friedr. Viehweg & Sohn, Braunschweig/ Wiesbaden, 
1989, in Unbehauen, "Regelungstechnik III" ["Control 
Engineering III"], 2nd edition, Friedr. Viehweg & Sohn, 

10 Braunschweig/Wiesbaden, 1986, in E. Pfeiffer, "Einfuhrung in 
die Dynamik" ["Introduction to Dynamics"], B.G. Teubner, 
Stuttgart, 1989, and in E. Ziegler (ed.), "Teubner-Taschenbuch 
der Mathematik" ["Teubner Manual of Mathematics"] B.G. 
Teubner, Stuttgart, 1996. When simulating, the problem 

15 frequently arises that a method for numerically solving a 

system of differential equations will terminate because the 
fundamental system of differential equations is singular. For 
the same reason, it can happen that an analog computer cannot 
foresee the system behavior. In addition, it frequently 

2 0 happens that solutions that are not plausible occur when 
solving the systems of differential equations. 

In such a case, it is known in the prior art that the 
programmer or a user of a program for solving differential 
25 equations is responsible for locating errors. With extensive 
experience, it is possible to locate the causes of the 
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singularity and, if appropriate, to reduce errors in the 
modeling of the technical system. This is generally very time 
consuming and costly. 

5 A method for automatic structure analysis of very large systems 
of differential-algebraic equations is described in S. Chowdry 
et al : "Automatic structure analysis of large scale 
differential algebraic systems' 7 IMTC 2001. Proceedings of the 
18th IEEE Instrumentation and Measurement Technology 

10 Conference, Budapest, Hungary, May 21-23, 2001, IEEE 

Instrumentation and Measurement Technology Conference (IMTC) : 
New York, NY: IEEE, US, Vol. 1 of 3 , Conf . 18, May 21, 2001 
(2001-05-21), pages 798-803. This method takes account both of 
numerical and of symbolic information, and also incorporates 

15 the effects of linear terms rigorously into the mathematical 
system, and is capable of providing a very accurate 
representation of the structure of a system of differential- 
algebraic equat ions . 

2 0 A structural algorithm that takes account of the structural 
properties of systems of differential -algebraic equations is 
described in the scientific article by J. Unger et al : 
"Structural analysis of differential-algebraic equation 
systems - theory and applications" Comput . Chem. Eng.; 

25 Computers & Chemical Engineering Aug 1995, Pergamon Press 

Inc., Tarrytown, NY, USA, Vol. 19, No. 8, August 1995 (1995- 
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08), pages 867-882. In that case, two systematically different 
structural approaches are implemented on a computer system in 
order to analyze systems of differential -algebraic equations 
structurally. These computer programs are applied both as 
5 separate programs for supporting the development of simulation 
models with a low index, and as constituents of the DIVA 
simulation environment, in order to facilitate a selection of 
suitable numerical algorithms, and to calculate consistent 
initial conditions . 

10 

The fact that the structural index of a linear differential - 
algebraic equation with constant coefficients and with the index 
1 can be arbitrarily high and to that extent contradicts a 
result previously published in the literature is shown by the 

15 scientific publication by Gunter Reissig et al : "Dif f erential- 
algebraic equations of index 1 may have an arbitrarily high 
structural index" Siam J. SCI. Comput.; Siam Journal of 
Scientific Computing 2000 SOC for Industrial & Applied 
Mathematics Publ . , Philadelphia, PA, USA, Vol. 21, No. 6, 2000, 

2 0 pages 1987-1990. This publication further states that when 

applied to differential -algebraic equations with an index of 1, 
Pantel ides' algorithm can execute an arbitrarily high number of 
iterations and differentiations. 
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Summary of the Invention : 

It is accordingly an object of the invention to provide a 
method for providing error information relating to 
inconsistencies in a system of differential equations which 
5 overcomes the above-mentioned disadvantages of the heretofore- 
known devices and methods of this general type. 

With the foregoing and other objects in view there is 
provided, in accordance with the invention, a method for 
10 predicting the behavior of a technical system starting from 
prescribed system properties, boundary conditions, and/or 
starting from prescribed influences on the system. The method 
comprises the following steps: 

providing a computer system ; 

15 defining a system of equations describing the technical 

system, the equations having the form f{{,x(t),x(t),... P x (k) (t),p)=0 , 
and including: 

20 : 

fn{t^t)^0.-^ k> (t),p}=0, 
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wherein x(t) and derivatives x(t),... t x (k) thereof in each case 
have m elements, and £ is a parameter vector that can occur in 
the system of equations; and 

executing a test method for providing error information 
5 relating to inconsistencies in the system of equations on the 
computer system, the test method having the following steps 1 
to 3 : 

step 1: setting up a dependence matrix A with m columns and n 
rows, and setting an element A to 

10 A(i,j)*0 when an i th row of f defined with f i (t,x(t) f x(t),...,x (k) (t)) is 
a function of 

a) a j th element of x expressed as Xj(t); or 

b) one of the derivatives of the j th element of x defined 
as Xj (s) (t) ; 

15 and otherwise setting the element A(i, j ) = 0 ; 

step 2 : determining a set of row ranks each having the numbers 
of those rows of the dependence matrix A that are mutually 
dependent and determining a set of column ranks of the 
dependence matrix A. each having the numbers of those columns 
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of the dependence matrix A that are mutually dependent if such 
row ranks and/or column ranks are present; and 
step 3: outputting error information containing the numbers 
contained in each row rank determined in step 2 and in each 
column rank determined in step 2 . 

The invention is based on the idea of examining the system of 
differential equations for inconsistencies with the aid of 
combinatorial methods as early as before a simulation, that is 
to say before executing numerical calculating steps or 
building an analog computer. When information relating to such 
inconsistencies is at hand, it is possible to deduce errors in 
the modeling of the technical system, or errors in the 
technical system itself, that can cause undesired 
contamination in a later simulation step, or even render the 
system incapable of being simulated. There is no need for a 
time-consuming execution of the simulation as far as 
termination, nor for a trial simulation. The correction of the 
model of the technical system or of the technical system 
itself is substantially facilitated by the presence of 
information relating to inconsistencies. 

The first step, of the method according to the invention, 
provides for setting up a dependence matrix A. The dependence 
matrix A has exactly as many columns as the dimension of the 
solution vector x of the system of differential equations. The 
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dependence matrix A has exactly as many rows as the number of 
differential equations in the given system of differential 
equations . 

5 The invention is not limited in this case only to systems of 
differential equations of the form f{t,x(t),x(t),.„,x (k) (t),p)=0. It 

can also be applied to special cases in which the system of 
equations for simulating the technical system has the form, 
for example, of f{t>x(t)) = 0 . In such a case, the talk is not of 

10 a system of differential equations but of a generally 
nonlinear system of equations. 

Moreover, it is not absolutely necessary for the system of 
differential equations that is provided for the simulation to 

15 be examined itself for inconsistencies. Moreover, another 
system of differential equations that describes the same 
technical system or the same technical process can be examined 
for inconsistencies. This can be advantageous when the last- 
named system of differential equations has a simpler structure 

2 0 than that to be simulated later. For example, in the 

simulation of electric networks the equations of modified node 
voltage analysis are usually the basis of the numerical 
simulation, although the so-called branch voltage-branch 
current equations have a substantially simpler structure. It 

25 is even possible frequently to determine more inconsistencies 
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or to limit inconsistencies more severely by examining systems 
of differential equations with a simple structure. 

The dependence matrix A in this case has elements A(i t j) that 

are set either to the value "zero" or to the value "nonzero" , 
without regard to the absolute values of the elements of A. In 
this case, the element A(i,j) can be set to a value * 0 when the 

i th row of f is a function of the j th element of X/ or of one 
of the derivatives of the j th element of x- l n a ll other 
cases, the element A(i,j) is set to the value "0" . 

It is immaterial for the method according to the invention 
whether the dependence matrix A is actually set up explicitly. 
Rather, what is important is the corresponding information as 
15 to whether, the i th row of £ is a function of the i th element of 
x or of one of its derivatives, or not, being available in the 
following method steps. 

In the second step of the method according to the invention, a 
2 0 set of row ranks is determined that in each case have the 

numbers of those rows of the dependence matrix A or those rows 
of the system of differential equations that are mutually 
dependent. Moreover, a set of column ranks of the dependence 
matrix A is determined that in each case have the numbers of 
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those columns of the dependence matrix A or those components 
of x or one of its derivatives that are mutually dependent. 

Here, "row rank" is understood as follows. A set C z of natural 
5 numbers i, /</<«, is a row rank of a matrix with n rows and m 
columns when it fulfils the following conditions: 

(i) there is no transversal T of the matrix A such that C z 
is contained in the set of row indices of T. 

10 

(ii) for each element c of C z there is a transversal T of A 
such that C z \{c} is wholly contained in the set of the 
row indices of T. 

15 Here, the expression C 2 \{c} represents that set which results 
when the element c is removed from the set C z . 

A transversal T of the matrix A is understood as follows. A 
transversal of a matrix A with n rows and m columns is one of, 
20 if appropriate, a plurality of possible sets of positions (i, 

j) of non-vanishing matrix entries A (i, j) , of which no two or 
more are in the same row or column. A set T of pairs (i, j) , 
l<i<n,l< j <m, constitutes a transversal of the matrix A if T 
fulfils the following conditions: 

25 
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(i) 



A (i, j) does not vanish for all elements (i, j) of T, 



and 



(ii) 



given that (i, j) and (i', j ') are two different 



elements of T, that is to say or j*f, it holds 



5 



that i*V and at the same time j * f . 



It is of less importance for the invention in this case 
whether one of the possible transversals of the dependence 
matrix A is calculated. The definition of the term 
10 "transversal" is required here, in order to illustrate the 
terms "row rank" and "column rank" . 

The term "row index of the transversal T" is understood as 
follows. The set Z of the "row indices of the transversal T" 
15 has as elements the row indices of the elements of T. In other 
words, this means that the set Z contains precisely those 
numbers i for which there is a j such that the element (i, j) 
is an element of T. 

2 0 A column rank is understood as a set C s of natural numbers i, 

l<i<m, of a matrix A with n rows and m columns when it fulfils 
the following conditions: 

(i) there is no transversal T of the matrix A such that C s 

2 5 is contained in the set of column indices of T. 
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(ii) for each element c of C s there is a transversal T of 
the matrix A such that C s \{c} is wholly contained in 
the set of the column indices of T. 

The u set of column indices of the transversal T" is understood 
as a set Z that contains as elements the column indices of the 
elements of the transversal T of the matrix A. In other words, 
this means that the set Z contains precisely those numbers j 
for which there is an i such that the element (i, j) is an 
element of the transversal T of the matrix A. 

After step 2 according to the invention is carried out, there 
is present a set of row ranks that respectively have the 
numbers of those rows of the dependence matrix A or those rows 
of the system of differential equations that are mutually 
dependent. Moreover, a set of column ranks of the dependence 
matrix A is present that in each case have the numbers of 
those columns of the dependence matrix £ or of those 
components of x or one of its derivatives that are mutually 
dependent . 

It is particularly easy having this information to deduce 
structural errors of the system of differential equations that 
is present as starting point. Such structural errors or 
inconsistencies in the original system of differential 
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equations are frequently causes of errors in the calculation 
of the solution of the system of differential equations. 



A basic step in executing the method according to the invention 
5 consists in finding row ranks and column ranks. Methods for 
finding row ranks and column ranks are known. Methods for. 
determining "row ranks" are to be found in the literature under 
the keyword of "minimally structurally singular subsets of 
equations" (compare C.C. Pantelides, The consistent 

10 initialization of differential -algebraic systems in SIAM J. 

Sci. Statist. Comput . , 9(2): 213-231, March 1988). In order to 
determine "column ranks", it is possible to use a method for 
determining "row ranks" when a transposed matrix is determined 
from that matrix of which the "column ranks" are to be 

15 determined. The row ranks of the corresponding transposed 

matrix are then determined in order to determine the column 
ranks of the matrix. 

This manner of determining the row ranks and column ranks for 
2 0 the method according to the invention is not to be understood 
in a limiting way. Rather, it is also possible to use other 
methods to determine row ranks and column ranks if the 
remaining conditions for the presence of row ranks and column 
ranks are fulfilled. 
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In the concluding step of the method according to the 
invention, for each row rank and for each column rank that had 
been determined in accordance with step 2, the numbers that 
are contained therein are output. In the case of the row 
5 ranks, these numbers indicate the running numbers of the 
equations of the system of differential equations that are 
possibly affected by a structural problem. In the case of 
column ranks, the numbers contained therein and output as 
error information indicate the numbers of the components of 
10 the solution vector x that are possibly affected by a 
structural problem. 

The method according to the invention can be used even before 
the execution of numerical calculating steps or the building 
of an analog computer to check the structure of the system of 
differential equations that is to be calculated for its 
consistency. If the method according to the invention had 
detected a system of differential equations as being 
inconsistent, it is possible to locate a large number of the 
possible structural errors of the system of differential 
equations. This accelerates the location of errors. In 
addition, time-consuming simulation tests are avoided. 

The method according to the invention can be applied for any 
25 type of a simulation system for the numerical solution of 
systems of differential equations, and for the solution of 
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systems of differential equations by means of an analog 
computer . 

In a particularly advantageous development of the invention, 
5 before executing step 1 of the method according to the 

invention an equation significance list of length n is applied 
in which each equation of the system of equations is assigned 
an equation number and/or an item of equation text 
information. In exactly the same fashion, before executing 

10 step 1 of the method according to the invention, a component 
significance list of length m is applied in which each 
component of the solution vector x is assigned a component 
number and/or an item of component text information. In this 
case, the equation text information or component text 

15 information respectively stored in the significance list is 

advantageously selected such that this information acquires a 
significance in connection to the technical system to be 
simulated. Consequently, constituents of the simulated 
technical system and substructures of the technical system to 

20 be simulated can be assigned equations of the system of 

differential equations and components of the solution vector x 
that facilitate an interpretation of the error information 
output with the aid of the method according to the invention. 

25 In step 3 of the method according to the invention, for this 
purpose it is provided to output the equation number and/or 
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the item of equation text information in accordance with the 
equation significance list instead of outputting the numbers 
contained in each row rank. For the case in which a whole 
number i is contained in the row rank, it is not the number i 
5 that is output but rather the contents of the i th component of 
the equation significance list. In exactly the same fashion, 
in step 3 of the method according to the invention the 
component number and/or the item of component text information 
in accordance with the component significance list is/are 
10 output instead of outputting the numbers contained in each 
column rank. For the case in which the whole number j is 
present in the column rank as a number, it is not the number j 
that is output, but rather the contents of the j th component 
of the component significance list. 

15 

Owing to this refinement of the method according to the 
invention, significance contents that refer directly to the 
technical system to be simulated are output as error 
information. This error information can be used to illustrate 
2 0 with particular ease, systematic errors in the structure of 
the system of differential equations for describing the 
technical system, such that locating the errors is further 
accelerated. 
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In the method according to the invention a preferably digital 
computer is used that has at least a memory, an arithmetic 
unit, an input device, and an output device. 

The invention is also implemented in a computer program for 
providing error information relating to inconsistencies in a 
system of equations. The computer program is designed in this 
case such that after inputting of the system properties, the 
initial or boundary conditions and the influences on the 
system, it is possible to execute an inventive method in 
accordance with one of the preceding claims. It is possible in 
this case to output results of a simulation, the solution 
vector or the solution vectors at different points in time. 
However, it is also possible for the computer program to 
provide only information relating to inconsistencies in 
accordance with one of the preceding claims. 

Because numerous defective program runs and/or simulation 
trials can be avoided, the computer program according to the 
20 invention results in substantial improvements in run time by 
comparison with the known simulation programs and simulation 
methods . 

Moreover, the invention relates to a data carrier with such a 
2 5 computer program, and to a method in which such a computer 

program is downloaded from an electronic data network, such as 
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from the Internet onto a computer connected to the data 
network . 

The computer system according to the invention is designed 
5 such that it is possible to execute on it a method according 
to the invention for providing error information relating to 
inconsistencies in a system of equations. 

Finally, the invention also relates to the use of a method 
10 according to the invention and/or a computer system for 

providing error information relating to inconsistencies in a 
system of equations. 

In the method according to the invention, the first step is to 
set up the system or system of differential equations that 
describes the process, or other forms such as, for example, 
nonlinear systems of equations that have no differential 
equations. It is assumed here, in this case, that the person 
skilled in the art is conversant with the mode of procedure 
required for this purpose. This is examined in detail in the 
specialist literature mentioned in the introduction to the 
description. There are also relevant computer programs for 
this purpose. 

25 Other features which are considered as characteristic for the 
invention are set forth in the appended claims. 
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Although the invention is illustrated and described herein, in 
general, as embodied in a method for providing error 
information relating to inconsistencies in a system of 
5 differential equations, it is nevertheless not intended to be 
limited to the details shown, since various modifications and 
structural changes may be made therein without departing from 
the spirit of the invention and within the scope and range of 
equivalents of the claims. 

10 

The construction and method of operation of the invention, 
however, together with additional objects and advantages 
thereof will be best understood from the following description 
of specific embodiments when read in connection with the 
15 accompanying drawings. 

Brief Description of the Drawings : 

Fig. 1 is a set of functions and a parameter vector 
illustrating a first exemplary embodiment of the invention; 

20 

Fig. 2 is a matrix and a corresponding set of transversal 
elements ; 

Fig. 3 is a set of functions and a parameter vector describing 
25 a second exemplary system; 
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Fig. 4 is a matrix, a transversal, and column and row results; 

Fig. 5 is an exemplary display of error results; 

5 Fig. 6 is a schematic of an equivalent circuit; 

Fig. 7 shows an equation significance list G and a component 
significance list K; and 

10 Fig. 8 is a resultant error list. 

Description of the Preferred Embodiments : 
The first exemplary embodiment in accordance with the 
invention, relates to the simulation of a technical system 
15 that can be described by the numerical solution of the system 
of equations 



20 



f(t,x(t),x(t),...,x (k) (t),p)=0 



of the form 



fl (*1 (0, *2 0)> * 3 (0, *1 (0, *2 (0. *3 (0, A > P 2 > />3 ) = 0 > 
fl Ol (0> *2 (0, *3 (0, *1 (0, *2 (0, *3 (0, Pi >Pl > Pi ) = o> 
A Ol (0. X 2 (0. *3 (0, *1 (0, *2 (0. *3 (0» Pi , P 2 , Pz ) = o, 



having the functions fl, f2 and f3, and having the parameter 
vector £ in accordance with Fig. 1. 
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In order to predict the behavior of the system, the system of 
equations is solved numerically, that is to say values for the 
unknown vector x(t) are calculated at one or more points in 
5 time t. For this purpose, known methods of numerical solution 
are used which run as a computer program on a computer system. 

According to the invention, steps for providing error 
information relating to inconsistencies in the system of 
10 equations are executed before the actual solution of the 
system of equations. 

The dependence matrix A specified in Fig. 2 is determined in 
step 1 of the method according to the invention, all non- 
15 vanishing elements of the dependence matrix A being denoted in 
Fig. 2 by a star ("*") . 

The element A(l,l) picked out by way of example is set to an 
essentially arbitrary non-vanishing value w *" , because the 

first element of f , that is to say f/t^t),gt),...,x ( * ) (t)), is a 

20 function of the first element of x, that is to say of Xi(t)- 

The element A(l,3) picked out by way of example is set to the 
value "0", because the third element of f, that is to say 
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f 3 (t,x(t),x(t),...,x (k) (t)), is independent of the first element of x, 

that is to say of x x (t) and of the derivatives of the first 
element of x# that is to say Xi (s) (t) . 

5 In step 2 of the method according to the invention, the first 

step is to determine the transversals T, specified in Fig. 2, of 
the dependence matrix A with the aim of determining row ranks 
and column ranks of the dependence matrix A specified in Fig. 2. 
Methods for determining transversals (of "maximum cardinality 

10 bipartite matchings") are specified in I.S. Duff, "On algorithms 
for obtaining a maximum transversal", ACM Trans. Math. Software, 
7 (3) :315-330, 1981, in E.L. Lawler, "Combinatorial Optimization: 
Networks and Matroids" , Holt, Rinehart and Winston, 1976, in L. 
Lovasz and M.D. Plummer, "Matching Theory", North-Holland 

15 Mathematics Studies 121. Annals of Discrete Mathematics, 29. 
North- Hoi land, 1986, or in C.C. Pantelides, "The consistent 
initialization of differential algebraic systems", SIAM J. Sci. 
Statist. Comput., 9 (2) : 213-231, Mar. 1988. 

20 The transversal T from Fig. 2 has 3 elements. Since n = m = 3, 
the dependence matrix A from Fig. 2 has neither row ranks nor 
column ranks . No such row ranks of the dependence matrix A are 
found by the method, specified in C.C. Pantelides "The consistent 
initialization of differential algebraic systems", for 

25 calculating row ranks ("minimally structurally singular 
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subsets of equations"). The method specified there for 
calculating row ranks also finds no such row rank of the 
transpose of the dependence matrix A. It follows that there is 
no column rank of the dependence matrix A. The set, determined 
5 in step 2 of row ranks of the dependence matrix A is empty. 
The set, determined in step 2 of column ranks of the 
dependence matrix A is likewise empty. Consequently, no error 
information is output in step 3 of the method according to the 
invention. 

10 

The system of equations specified at the beginning can be 
solved with a higher probability without the occurrence of 
errors. If the simulation is terminated nevertheless, or if it 
supplies implausible solutions, or if it is entirely 
15 impossible, then it is possible to exclude structural errors 

in the description of the technical system or in the technical 
system itself with a high degree of probability. This 
substantially simplifies the location of errors. 

2 0 The second exemplary embodiment in accordance with the 

invention relates to the simulation of a further technical 
system (not shown here) , whose behavior can be described by 
the numerical solution of the system of equations 



25 
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of the form 

fl (*1 (0> X 2 (0. *3 (0» A > P 2 > />3 ) = °> 
fl (*\ (0» X 2 (0, * 3 (0, Pl,P 2 >P3) = °> 
h (*l (0, x 2 (0. *3 (0, A > />2 > />3 ) = °> 

5 

having the functions fl, f2 and f3 and having the parameter 
vector £ in accordance with Fig. 3. 

It is aimed to solve the system of equations numerically in 
10 order to predict the behavior of the system, that is to say 

values are to be calculated for the unknown vector x(t) at one 
or more points in time t. Use is made for this purpose of 
known methods of numerical solution (not illustrated here) 
that run as a computer program on a computer system (not shown 
15 here) . 

According to the invention, methods for providing error 
information relating to inconsistencies in the system of 
equations are executed for the actual solution of the system 
20 of equations. 

The dependence matrix A specified in Fig. 4 is determined in 
step 1 of the method according to the invention, all non- 
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vanishing elements of the dependence matrix A being denoted by 
a star ("*") in Fig. 4. 



The element £(1,1) picked out by way of example is set to an 
5 essentially arbitrary non-vanishing value , because the 

first row of f, that is to say fi{f,x(t) 9 x(t) 9 ...,x {k) (t)) , is ^ function 

of the first element of x, that is to say of Xi(t) • The same 
holds for the elements A (1,2), A (1,3), A (2 ,3) and A ( 3 , 3 ) . 



10 The element A (3,1) picked out by way of example is set to the 
value u 0", the third element of £_, that is to say 
Aif»x{t) 9 x(t) 9 ... 9 x w (t)) § is independent of the first element of x, 

that is to say of Xi(t) and of derivatives of the first 
element of X/ that is to say Xi (s) (t) . The same holds for the 
15 elements A(2,l) , A(3,2) and A(2,2) . 



In step 2 of the method according to the invention, the first 
step is to determine the transversals T, specified in Fig. 4, of 
the dependence matrix A with the aim of determining row ranks 
2 0 and column ranks of the dependence matrix A specified in Fig. 4. 

The row rank {2, 3} of the dependence matrix A is found by the 
method for calculating row ranks ("minimally structured 
singular subsets of equations") specified in C. C. Pantelides 
25 "The consistent initialization of differential algebraic 
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systems''. The set Z of the row ranks, found in step 2, of the 
dependence matrix A is specified in Fig. 4. 

The row rank {l, 2} of the transpose of the dependence matrix 
A is found by the method, specified in C.C. Pantelides, "The 
consistent initialization of differential-algebraic systems", 
for calculating row ranks applied to the transpose of the 
dependence matrix A. According to the invention, this row rank 
of the transpose of the dependence matrix A is understood as 
the column rank of the dependence matrix A. The set S of the 
column ranks, found in step 2, of the dependence matrix A is 
specified in Fig. 4. The error information specified in Fig. 5 
is output in step 3 . 

In accordance with the invention, there is not even an attempt 
to predict the behavior of the basic system using numerical 
means, because in this case errors will occur if a simulation 
is at all possible. However, the modeling of the system and 
the system itself need to be checked once again. This saves 
valuable computing time on the computer system (not shown 
here) . Error locating is substantially simplified by the 
knowledge of the error information output in step 3. 

A third exemplary embodiment in accordance with the invention 
relates to the technical system shown in Fig. 6, whose 
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behavior can be described by the numerical solution of the 
system of equations 

/(r 3 x(0,2(0v.^ ( * ) (0,£)=0 / 

5 

of the form 

fx (*i (0, *2 (0, * 3 (0, Pi , Pi >P 3 ) = °> 
fi (*i (0, * 2 (0, * 3 (0> Pi s P 2 > />3 ) = °> 

/ 3 (*j (0, * 2 (/), * 3 (0, />! > />2 > />3 ) = °> 

10 having the functions fl, f2 and f3 and having the parameter 
vector £ in accordance with Fig. 3. 

It is aimed to solve the system of equations numerically in 
order to predict the behavior of the system, that is to say 
15 values are to be calculated for the unknown vector x(t) at one 
or more points in time t. Use is made for this purpose of 
known methods of numerical solution that run as a computer 
program on a computer system. 

2 0 Error information relating to inconsistencies in the relevant 
system of equations is provided in this case as follows 
according to the invention. 
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The solutions of the relevant system of equations are 
quiescent states or ''operating points" or U DC solutions" of 
the electric network specified in Fig. 6 and which comprises 
the following network elements ("components"): 

5 ° a linear resistor with resistance value R between the 

nodes 1 and 2 of the network, 

D a linear capacitor with a capacitance CI between the 
nodes 1 and 0 of the network, and 

° a linear capacitor with a capacitance C2 between the 
10 nodes 2 and 0 of the network. 



This results in the components of the parameter vector p, 
specified in Fig. 3, being according to the values CI, 1/R, 
C2 . 

15 

The components xl(t), x2(t) and x3 (t) of x( fc ) correspond to 
the following variables of the network specified in Fig. 6: 

° xl(t) corresponds to the voltage between the nodes 1 and 
0, 

2 0 0 x2 (t) corresponds to the voltage between the nodes 2 and 

0 , and 

a x3(t) corresponds to the voltage between the nodes 1 and 
2 . 
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The first equation of the system of equations, that is to say, 
the equation 

fl(xl(t), x2(t), x3(t), pi, p2, p3) = 0 

5 

is the Kirchhoff voltage equation for the mesh comprising all 
the three network elements of the network from Fig. 6. 

The second equation of the system of equations, that is to 
10 say, the equation 

f2(xl(t), x2(t), x3(t), pi, p2, p3) = 0 

is the Kirchhoff current equation for the nodes 1 of the 
15 circuit in Fig. 6. 

The third equation of the system of equations, that is to say, 
the equation 

20 f3(xl(t), x2(t), x3(t), pi, p2, p3 ) = 0 

is the Kirchhoff current equation for the nodes 2 of the 
circuit in Fig . 6 . 
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The equation significance list G specified in Fig. 7 and the 
component significance list K specified in Fig. 7 are designed 
in accordance with the invention. 

The set Z specified in Fig. 4, of the determined row ranks, 
and the set S, specified in Fig. 4 of the determined column 
ranks are determined as in the preceding exemplary 
embodiments . 

The error information specified in Fig. 8 is output in step 3 
of the method developed further in accordance with the 
invention, use being made of the equation significance list G 
and the component significance list K from Fig. 7. 

It may be seen from the error information specified in Fig. 8 
that the Kirchhoff current equations relating to the nodes 1 
and 2 of the network from Fig. 6 are linearly dependent on one 
another, and that the two voltages between nodes 1 and 0 and 
between nodes 2 and 0 of the network from Fig. 6 for quiescent 
states of this network cannot be uniquely determined. 

In accordance with the invention, there is not even an attempt 
to predict the behavior of the basic system using numerical 
means, because in this case errors will occur if a simulation 
is at all possible. However, the modeling of the system and 
the system itself need to be checked once again. This saves 
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valuable computing time on the computer system. Error locating 
is substantially simplified by the knowledge of the error 
information output in step 3 . 

5 The term "computer-readable medium" as used herein should be 
understood in its broadest sense. That is, it includes, at the 
least, any kind of computer memory such as floppy disks, hard 
disks, CD-ROMs, flash ROMs, non-volatile and volatile ROM and 
RAM, and memory cards, as well as carrier signals for distance 
10 communication. 



-31- 



